Dynamical dimensional reduction in toy models of 4D causal quantum gravity 
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In recent years several approaches to quantum gravity have found evidence for a scale dependent 
spectral dimension of space-time varying from four at large scales to two at small scales of order 
of the Planck length. The first evidence came from numerical results of four- dimensional causal 
dynamical triangulations (CDT) [Ambj0rn et al., Phys. Rev. Lett. 95 (2005) 171]. Since then 
little progress has been made in analytically understanding the numerical results coming from the 
CDT approach and showing that they remain valid when taking the continuum limit. In this letter 
we propose a new toy model of "radially reduced" four-dimensional CDT in which we can take the 
continuum limit analytically and obtain a scale dependent spectral dimension varying from four to 
two with scale. Furthermore, the functional behaviour of the spectral dimension is exactly of the 
form which was conjectured on the basis of the numerical results. 



Introduction - The quest to reconcile classical gen- 
eral relativity and quantum mechanics has a long history. 
One of the main difficulties is the fact that gravity is 
perturbatively non-renormalizable in four dimensions as 
was shown by 't Hooft and Veltman in the seventies [l|. 
However, defining quantum gravity non-perturbatively, 
there is evidence from different approaches that there 
might still be a non-trivial ultraviolet fixed-point as sug- 
gested by Weinberg [2| . The causal dynamical triangula- 
tion (CDT) approach to quantum gravity (see Q and Q 
for a review) recently gave a surprising answer to what 
dynamical mechanism might regulate the theory at short 
distances in such a scenario. In particular, numerical sim- 
ulations [E[ show evidence for the dynamical reduction of 
the so-called spectral dimension from four at large scales 
to two at small scales of order of the Planck length (see 
also [6] for a discussion in a three-dimensional setting). 
More recently similar results have also been observed in 
other approaches to quantum gravity, most notably the 
exact renormalization group approach [7] and Horava- 
Lifshitz gravity @ . This also points towards several pos- 
sible relations between those approaches [|[. 

Previous numerical insights from four- 
dimensional CDT - The idea behind the spectral 
dimension as a scale dependent measure of dimension- 
ality is the following. Consider a diffusion process on a 
fixed space-time geometry. The diffusion kernel K g is 
determined by the diffusion equation 
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K g (y, y ; cr) = A g K g (y, y ; cr), 



(1) 



where g denotes the space-time metric, yo the starting 
point of the diffusion and y the position of the diffusion 
process after time a. We can then define the return prob- 
ability (density) by choosing starting and end point equal 
and integrating over all positions 



p g(°) 
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d d ys/detg ab (y) K g (y,y;a) (2) 



with V g = J d d y y^det g a b(y)- If we consider diffusion on 
quantum space-time we would have to take the ensem- 
ble average, formally defined using the gravitational path 



integral 



| J V[g ab ] e- s »^P g (v), 



(3) 



where Z = / T>[g ab ] e 



-S E (g a b) 



is the partition function 



and Ssigab) the Euclidean Einstein Hilbert action. 

Formallyjthe sale dependent spectral dimension is then 
defined as [5j 



D s (a) 



dlog(P(a)) z 



(4) 



d log a 

where a is the diffusion time and corresponds to the scale 
at which the diffusion process probes the quantum geom- 
etry. 

In the CDT approach the gravitational partition func- 
tion is defined as a sum over four-dimensional simplicial 
geometries with a strict time-sliced structure 



S E (T) 



(5) 



where Ct is a combinatorial symmetry factor and Se(T) 
the Euclidean Einstein- Regge action of the triangulation 
T; S E (T) = XN 4 (T) - vN 2 (T). Here A is the bare cos- 
mological constant and v is the bare inverse Newton's 
constant. N^T) denotes the number of four-simplices 
and N 2 (T) the number of two-simplices in T. 

In theory, one would like to evaluate the partition func- 
tion (|5]) and take the infinite volume limit in which A 
is tuned towards its critical value and v is expressed in 
terms of the inverse renormalized Newton's constant 1 /G. 
While an analytical evaluation of (j3j) in four-dimensional 
CDT, i.e. using the partition function ([5]), is out of reach, 
Monte Carlo simulations of random walks (discrete dif- 
fusion) on four-dimensional CDT [E| yield a scale depen- 
dent spectral dimension given by 



D s (a) =4.02 



119 I 1.80 ±0.25, cr^O 
54 + a I 4.02 ±0.1, cr^oo. 



(6) 



A possible problem with (j6j) is that the simulations are 
inevitably affected by finite size effects and the dimen- 
sional reduction observed might simply be an artefact of 




Figure 1. The two building blocks of four-dimensional CDT; 
The (4,1) and (3,2) simplex are on the left and right respec- 
tively. Dotted lines correspond to 3 simplices. 



the discreteness scale. However, if one assumes that this 
expression can be extrapolated to continuum physics, one 
would conjecture the following expression for the return 
probability (density) © for four-dimensional CDT in the 

infinite volume limit 
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a 2 1 + const. G/a' 



(7) 



In this letter we introduce a new toy model of "radially 
reduced" four-dimensional CDT, a particular so-called 
multigraph ensemble in which the spectral dimension can 
be calculated analytically. At this point we would like 
to remind the reader of the difficulty caused by the two 
kinds of randomness - we are describing a random walk 
on a random geometry. The idea behind the "radial re- 
duction" is that what essentially determines the prop- 
erties of the random walk are the volume growth and 
the resistance growth (viewing the graph as an electri- 
cal network [lQj)- Hence, by omitting degrees of freedom 
which do not influence these quantities we get an effec- 
tive description of the random walk on the original CDT 
ensemble. 

Many of the more technical aspects of the present work, 
such as mathematically rigorous proofs, will be left for a 
forthcoming publication (ill ] . 

Multigraphs as toy models for causal quantum 
gravity - Consider four-dimensional rooted causal tri- 
angulations [l2| of topology I x S 3 which are made of 
(4,1) and (3,2)-simplices connecting vertices at distance 
n to vertices at distance n + 1 from the root (Figure [1]). 
Let us denote by T(t) the ball of radius t around the 
root. This subset of the triangulation is characterised 
by its bulk variables N^t) = Ni(T(t)), i = 0,1,2,3,4 
which denote the number of i-simplices in this section of 
the triangulation. We can further distinguish these vari- 
ables; for example, there are two different four-simplices 



so7V 4 (£) =N^ 1 \t)- 
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AQ°' ; (t). Further, there are three 



different types of three-simplices and two different types 



1 Note that throughout this letter we write 
multiplicative logarithmic corrections, while 
ity up to a multiplicative constant. 
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Figure 2. An example of a multigraph. 



of triangles and links (i.e. space-like and time-like). These 
ten different bulk variables are related by seven topolog- 
ical constraints, leading to only three independent vari- 
ables 0. 

We can now define a multigraph from a given causal 
triangulation by the injective map in which we shrink all 
vertices at a fixed graph distance n from the root to a 
single point (see Figure [2]). Let us denote the number 
of time-like edges connecting vertex n to vertex n + 1 
by L n . The multigraph M is then completely charac- 
terised by the sequence M = {L ni n = 0,1,2,...}. Fur- 
thermore, the partition function (|5]) induces a measure 
u(v) for the {L ni n = 0, 1, 2, ...} which defines the multi- 
graph ensemble. In the infinite volume limit this measure 
n{y) depends only on the bare Newton's constant l/v as 
an external parameter. 

In previous work on two-dimensional CDT [HI it was 
observed that multigraph ensembles provide us with a 
"radial" approximation of the causal triangulation. In 
particular, it appears that studying random walks on 
these ensembles captures the full behaviour of the spec- 
tral dimension in two-dimensional models of CDT [13| . 
This suggests that they might also serve as a realis- 
tic toy model for studying the spectral dimension of 
higher-dimensional CDT. However, the work [l3| on two- 
dimensional CDT does not describe a scale dependent 
spectral dimension. In a forthcoming article [ll| we show 
that by introducing a length scale in the model one can 
actually obtain a spectral dimension varying from two at 
large scales to one at small scales. This extends a toy 
model of scale dependent spectral dimension using ran- 
dom combs first described in [l4| and can be extended to 
four-dimensional CDT, as we describe here. 

While an analytical derivation of the measure u(y) 
of the multigraph ensemble related to four-dimensional 
CDT is still absent, we know for example from computer 
simulations [15[ that (N^t)) z ~ t 4 . From the topologi- 
cal relations mentioned above we find for the number of 
time-like links N± L (t) (up to boundary terms) that 



N- 



TL 



2N (t) + ±N ( *' 2) (t) 



(8) 



and using N (t) < N 4 (t)/5 and N^ 2 \t) < N 4 (t) we get 



(Ni(t)) z < const. (N 4 (t)), 



(9) 



This suggests that for the measure u(v) of the multi- 
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graph ensemble related to four-dimensional CDT 

< b (*)U) = (EM ( 10 ) 

\n=0 I 

Motivated by an analytical model of a multigraph en- 
semble with scale dependent spectral dimension derived 
from two-dimensional CDT, as mentioned above [ll|, we 
are led to the following ansatz on the expectation of the 
number of links 



N , 



vN 



■N 



(11) 



where e > can be taken arbitrarily small. The use of e 
here is for purely technical reasons and for all practical 
purposes one can think of it as being zero. This leads to 
the volume expectation of the form 

cH M <(B(% M <cM pM (12) 

where c, c are positive constants. Further, we assume that 
the resistance from a vertex at distance N to infinity, i.e. 
R(N) = Y^=n r~> anu - connectivity at height N, i.e. 
Ljv, are bounded above by 



R(N) < 



N 



Ln < (Ln)^ HV^N 



) 



(13) 
(14) 



for N > No for almost all graphs of the ensemble (such 
as those obtained from self- averaging in the computer 
simulations) where ip(x) and ijj+(x) are functions which 
diverge and vary slowly at x = and x = oo. For 
example in the two dimensional model we know that 
ip(x) ~ ^-\-(x) ~ I log x I [l3|| . These rather technical as- 
sumptions (fT3|) and JT3J) which allow only slowly varying 
fluctuations around the average are actually quite nat- 
ural as one can see in the lower- dimensional study [TT| . 
Note that the finite resistance R(N) implies that this 
multigraph ensemble is non-recurrent, i.e. d s > 2. 

Spectral dimension and dimensional reduc- 
tion - We propose the multigraph model with measure 
u(v) defined by (fTT]h (fT3|) and (fH|) as a toy model for 
four-dimensional CDT. In this section we show that it 
has a scale dependent spectral dimension which is two at 
small scales and four at large scales. 

Let pM(t) be the probability that a random walker on 
a fixed multigraph M returns to the root after time t. 
The spectral dimension is defined through pm{P) y i a the 
relation pm($) ~ t~ ds l 2 at large time. More precisely, 
using the generating function for the return probabilities 



Qm(x) y^PM^)(l - X) 



t/2 



t=0 



we define the spectral dimension using 

Qm(x) ~ x- 1+ds/2 as x^O 



(15) 



(16) 



in the case where the random walk is recurrent (d s < 2) 
and Qm(x) diverges as x — >> 0. If the random walk is 



non-recurrent, Qm(x) is finite and we define the spec- 
tral dimension through the derivative of Qm(x) of lowest 
degree which is diverging via the relation 



x -)• 



(17) 



for 2k < d s < 2(k + 1). As we will see later for u(y) 
the expectation value of Qm(x) is finite while its first 
derivative diverges. 

Given a fixed multigraph M the probability for a ran- 
dom walker at n to step next to n + 1 is given by 
p n (M) = L n /(L n -i + L n ) and the probability that the 
next step is to n — 1 is 1 — p n (M) (note that the proba- 
bility to move from the root to vertex one is 1). Then we 
decompose the random walk into two pieces; a step from 
vertex nton + 1, then a random walk returning to n + 1 
and a final step from n + 1 to n at time t. This decompo- 
sition relates QM n (x) and QM n+1 (x) and the generating 
function satisfies the following recursion relation [13| 

VM n (x) rjM n+1 (x) + y xL n f] Mn (x)r]M n+1 (x) (18) 

where r]M n (x) = QM n (x) / L n and M n is the multigraph 
obtained from M by removing the first n vertices and 
all edges attached to them and relabelling the remaining 
multigraph. Recall that Qm(x) = Qm ( x ). Differentiat- 
ing (fT8|) and iterating we get 



N-l 



Wm (x)\ = Wm n (x)\-1[ ^ 



xL n f] Mri (x) 



n=0 



" xL n 7] Mn+1 (x) 



(19) 
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L nVM n (x)VM n+1 ( X ) 
+ xL n f] Mn+1 (x) 



n—1 1 



xL k f] Mk (x) 



which is the starting point of our proofs. Starting from 
CEJ) and using (fT3]) - (p^|) we find that 



'■x^BtN^ < Wm ° {X)1 < 



/ [ VM n * (0) 



N* 



^ L n VM n (0)VM n+1 (0) 



(20) 



n=0 



where N* = [for -1 / 2 ] and c, d and b are positive 
constants. The detailed derivation of (1201) is more in- 
volved and will be presented elsewhere [11]. Noting that 
Vm n (0) = R(N) is the (finite) resistance, averaging over 
the ensemble u(y) and using (fTT]) - (fT4|) and Jensen's in- 
equality we have that 



1 



12 v - 



\V^X 2 + 4] 

< (IQ , m(^)IU)< c + -4 - -(21) 



V + X 



We obtain the scaling limit by taking the lattice spac- 
ing a — >• and following the prescription of [14] define 

\Q'{d,G)\ ee lim (fL) (\Q' M (x = Ol)^) (22) 



with v 



1 



2 /G. Using ([20]) we now have that 



<|Q'(C,G)|<c + 



j>l{G-i^-i+i) 



.(23) 



Eq. ([23]) implies that <i s = 2 in the short walk limit 
(i.e. £ — >• oo ) and d s = 4 — e in the long walk limit (i.e. 
£ — >• 0) for e arbitrarily small (see for the detailed cor- 
respondence between walk length and £). At this point 
we should mention that it is a highly non-trivial fact that 
there exists a non-trivial limit ([22]) . The first example of 
the existence of such a limit (in the recurrent case) was 
given in the context of random combs [14j . 

We can extract from ([2T]) the average return proba- 
bility as a function of large walk length. In particular, 
using a Tauberian theorem one gets (setting e to zero in 
the expressions below to simplify the discussion) 
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l 



(24) 



as t — y oo. Scaling t(a) — [cr/aj and v(a) = a/G as be- 
fore one obtains the probability density of the continuous 
diffusion time a through 



>) = 1™o(§)~ 1(pmW W) 



2G 2 



a 2 1 + 2G/ct 



.(25) 



This is precisely expression (0 which was conjectured in 
[E[ as the behaviour of the continuum return probability 
density for diffusion on four-dimensional CDT. It yields 
the scale dependent continuum spectral dimension 



D s (a)=4 1 



1 



2 + a/G 

consistent with the numerical results. 



(26) 



Discussion - We present a multigraph ensem- 
ble which describes "radially reduced" four-dimensional 
CDT. It is shown that starting from ([Til), $M and (HS> > 
which are all motivated from lower-dimensional studies 
one can calculate the spectral dimension of the re- 
sulting model analytically (see also for many of the 
technical details of the derivation). One observes a scale 
dependent spectral dimension varying from four at large 
scales to two at small scales. In fact, one of the non- trivial 
results of this letter is that one can actually perform the 
continuum limit exactly and obtain an expression for the 
continuum return probability density ([25]) . This result is 
compatible with the numerical results from [5| and fur- 
thermore gives a derivation of the conjectured expression 
([7]) . It is quite surprising that the above assumptions suf- 
fice to derive this result. While ([TT]) only describes the 
expectation value of time-like edges at distance N from 
the root it would be interested to further study correla- 
tions of the L n from the numerical data and possible de- 
rive their exact distribution. In particular, such a study 
should shed further light on the mechanism of dynami- 
cal dimensional reduction in four-dimensional CDT and 
would be a first step towards an analytical solution of the 
full four-dimensional model. 

Acknowledgements - GG would like to acknowl- 
edge the support of the A.G. Le vent is Foundation and 
the A.S. Onassis Public Benefit Foundation grant F-ZG 
097/ 2010-2011. This work is supported by EPSRC 
grant EP/I01263X/1 and STFC grant ST/G000492/1. 
SZ would like to acknowledge support of the STFC, as 
well as of a EPRSC visiting grant and Visiting Scholar- 
ship at Corpus Christi College, Oxford University during 
which this work was started. 



[1] G. 't Hooft and M. J. G. Veltman, Phys. Theor. A20 
(1974) 69-94. 

[2] S. Weinberg, in General Relativity, S. Hawking and W. 

Israel, eds., p. 790-831, 1980. 
[3] Ambj0 rn, J., Loll, R. N ucl. Phys. B536, 407-434 (1998), 

arXiv: |hep-th/9805108| 
[4] J. Ambj0rn, A. Gorlich, J. Jurkiewicz, R. Loll, 

[arXiv: 1007.25601 [hep-th]] . 
[5] Ambj0rn, J., Jurkiewicz, J., Loll, R. Phys. Rev. Lett. 

95, 171301 (2005), arXiv: |hep-th/0505Tl3 
[6] D. Benedetti, J. H enson, Phys. Rev. D80 (2009) 124036. 

[arXiv:0911.040T1 [hep-th]] . 
[7] D. F. Litim, Ph ys. Rev. Lett. 92 (2004) 201301, arXiv: 

|hep-th/0312114| O Lauscher and M . Reuter, JHEP 10 

(2005) 050, arXiv: |hep-th/0508202l 
[8] P. Hofava, Phys. Rev. D79 (2009) 084008, arXiv: 

0901.3775 [hep-th]; Phys. Rev. Lett. 102 (2009) 161301, 

arXiv: 0901.3775 [hep-th]. 
[9] J. Ambj0rn, A. Gorlich, S. Jordan, J. Jurkiewicz, R. Loll, 

Phys. Lett. B690 (2010) 413-419. |arXiv:1002.3"298l [hep- 



th]]; J. Ambj0rn, S. Jordan, J. Jurkiewicz, R. Loll, 
[arXiv: 1108.39321 [hep-th]]; T. P . Sotiriou, M. Visser, 
S. Weinfurtner, [arXiv: 1 105.6098 [hep-th]]; C. Anderson, 
S. Carlip, J. H. Cooperman, P. Horava, R. Kommu and 
P. R. Zulkowski, "Quantizing Horava- Lifshitz Gravity 
via Causal Dynamical Triangulations," arXiv:1111.6634 
[hep-th]. 

[10] Lyons, R. with Peres, Y. (2011). Probability on Trees and 
Networks. Cambridge University Press. In preparation. 

[11] G. Giasemidis, J. F. Wheater and S. Zohren, (2012) to 
appear. 

[12] J. Ambj0rn, J. Jurkiewicz, R. L oll, Nucl. Phys. B610 

(2001) 347-382. [hep-th/0105267] . 
[13] B. Durhuus, T. Jonsson and J. F. Wheater, J. Stat. Phys. 

139 (2010) 859-881, larXiv : 0908. 3643 ] [mat h-ph]; J. St at. 

Phys. 128 (2007) 1237-1260, |arXiv:math-ph/0607020] 
[14] M. R. Atkin, G. Gia semidis, and J. F . Wheater, J.Phys. 

A44 265001 (2011). la7Xiv:1101.4174l [hep-th] 
[15] J. Ambj0rn, J. Jurkiewicz, and R. Loll, Phys. Rev. Lett. 

93 (2004) 131301, arXiv: hep-th/0404156| 



